Thermodynamic magnetization of a strongly interacting two-dimensional system 
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We report thermodynamic magnetization measurements of a 2-dimensional electron gas for several 
high mobility Si-MOSFETs. The low-temperature magnetization is shown to be strongly sub-linear 
function of the magnetic field. The susceptibility determined from the zero-field slope diverges 
as l/T™, with a = 2.4 ± 0.2 even at high electron densities, in apparent contradiction with the 
Fermi-liquid picture. 

PACS numbers: 



Magnetic properties of the strongly-interacting elec- 
tron gas have long been a subject of intensive theoretical 
and experimental investigations. The Coulomb interac- 
tion, as reflected by the exchange term, favors parallel 
spins, and therefore leads to ferromagnctism, whenever 
it is strong enough. The strength of the interaction is de- 
termined by the ratio between the typical Coulomb and 
kinetic energies, and is customarily characterized by a 
dimensionlcss parameter r s = p/ag, ag being the Bohr 
radius in the material, and p ~ n~ x l d is the effective dis- 
tance between electrons in d dimensions; in two dimen- 
sions p = l/y/irn, n being the electron density. When n 
decreases, the relative effect of the interactions becomes 
stronger, contrary to naive expectations. At zero temper- 
ature and sufficiently low density the clean single-valley 
system of itinerant electrons was predicted to become fer- 
romagnetic, a phenomenon called Stoner instability [l[ . 
In two dimensions, according to the Mcrmin- Wagner the- 
orem, finite temperature destroys the ferromagnetic or- 
der. 

In a real system there is always some degree of disor- 
der, which favors an antifcrromagnctic interaction. In- 
deed, the ground state of two localized spins is a singlet, 
much like a hydrogen molecule. A system of many lo- 
calized spins preserves the tendency to order neighbor- 
ing spins in opposite directions, that is, the coupling is 
antiferromagnetic; see 0] for a review. Intensive investi- 
gations of magnetic properties of doped semiconductors, 
particularly phosphorus doped Si, in the 1980's led to 
a substantial understanding of the interplay between in- 
teractions and disorder. The experimentally observed di- 
vergence of the susceptibility as l/T a , with a m 0.6, on 
the low-density side of the metal-insulator transition has 
been well-understood ; however the persistence of the 
divergence on the high-density side remains a puzzle. 

Two-dimensional gated structures like MOSFETs al- 
low for a gradual change in electron density, and hence 
interaction strength, without strongly affecting sample 
disorder. Interest in the magnetic properties of two- 
dimensional electron gas (2DEG) was stimulated by the 



observation of strong suppression of conductivity in Si 
MOSFETs by an in-plane magnetic field which po- 
larizes electrons and drives the system into an insulat- 
ing phase. Scalinganalysis of the magnetoresistance led 
the authors of0, TH to suggest a quantum phase tran- 
sition into a ferromagnetic state at the density of the 
metal-insulator transition, n c Q . This claim was con- 
tested in [?[ and [1] on the basis of thermodynamic and 
Shubnikov-de Haas measurements, respectively. 

In this Letter we use the technique developed and 
described in Q to study the magnetization at elevated 
temperatures. This technique measures the recharging 
of the gate-to-2DEG capacitor due to a change in the 
2DEG chemical potential pi. By modulating an external 
magnetic field with the amplitude 8B{uj), while keeping 
the gate voltage constant we are able to measure a gate 
charge modulation 8Q(uj) given by 
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C{uj) dp 
e dB 
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where the capacitance C(ui) is measured independently 
by modulating the gate voltage at the same frequency. 
We note that all the quantities in Eq. [1] also depend on 
the electron density n; this dependence is not explicitly 
indicated. By virtue of the Maxwell relation dp/dB = 
—dm/dn, dp/dB can be expressed as the derivative of 
the magnetization per unit area m with respect to the 
density. In principle, integrating dm/dn over the density 
from n = would give m(n). 

In this Letter we present measurements from deep into 
the insulating region to far into the metallic phase, al- 
lowing us to integrate dm/dn from almost n — 0. The 
magnetization m(n) thus obtained exhibits a surprisingly 
strong temperature dependence at low magnetic field 
pgB <SC T even at high densities. This result contra- 
dicts expectations from the Fermi liquid theory. Indeed, 
at high densities interactions are relatively weak, there- 
fore magnetization should be temperature independent at 
temperatures small compared to the Fermi energy Ep. 

Extension of the measurements deep into the insulat- 
ing regime became possible due the lower sample resis- 
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tance at elevated temperatures, together with good qual- 
ity of the contacts. It is important to note that the 
expression ([T]) holds even if the capacitance drops and 
acquires imaginary part due to the contact and chan- 
nel resistance, which happens at the resistance about 
500 MOhms. This facilitates measurements deep into 
insulating region down to less the half the density of the 
metal-insulator transition in the sample. 

Our measurements were performed on several high- 
mobility Hall bar shaped Si-MOS structures with 2DEG 
located on the (100) interface between Si surface and 
SiC>2. Such a 2DEG possesses a two-fold valley degen- 
eracy in addition to its spin degeneracy. We used two 
groups of samples: those fabricated in Russia (R), similar 
to those used in Ref.EI, 0, Q and those made in Holland 
(H) , used in Ref. [1, [J The in-plane magnetic field was 
modulated at frequency lo/Itx = 6.1 Hz with amplitude 
40 mT. Standard He4 pumping and heating were used 
to set the temperature in the range 1.7-13K. We present 
the results obtained from the R sample with peak mo- 
bility 3.4m 2 /Vs and H sample with 3.3m 2 /Vs at 1.7K, 
measured most extensively; the data for several other R 
and H samples were very similar. An example of the data 
collected at different temperatures and electron densities 
is shown in Fig.[T] At low field, dm/ dn(B) grows linearly 
below some temperature-dependent field B* indicated by 
the gray area in Fig.[TJ above which the slope changes. 
The maximal dm/dn{B) at lowest achievable density and 
temperature reaches 0.9/is at B ss 0.3 T, indicating al- 
most full spin polarization. 

Let us first consider the low-field slope of the data, 
d 2 m/dndB\B=o = dx/dn\B=o, plotted in Fig. [5] as a 
function of density at different temperatures. The slope 
is large and positive at low temperature and density. As 
the density increases, the slope decreases, changes sign, 
and finally almost tends to zero at highest densities. As 
a function of temperature the slope decreases rapidly al- 
most to zero at T — 13K. There is some small negative 
slope left at highest temperatures in the R samples, which 
we attribute to the diamagnetic contribution due to the 
finite thickness of the 2DEG Q ■ This negative slope is 
absent for H samples [l(| • 

In order to obtain the zero field susceptibility we inte- 
grate dx/dn: 

n 

X (n,T) = J d X /dn(n',T)dn' (2) 
o 

We were able to measure dx/dn down to densities < 
iil(T) < n c , so we must extrapolate the integrand in @ 
down to zero density. For simplicity we took dx/dn to 
be constant and equal to its value at nx,(T). In order 
to get the spin part of susceptibility we need to subtract 
the temperature-independent diamagnetic susceptibility. 
This we did by equating the high temperature suscepti- 
bility to the non-renormalizcd Pauli one. The results of 




in-plane magnetic field B (Tesla) 

FIG. 1: Differential magnetization, dm/dn, for several den- 
sities n are given in 10 n cm~ 2 . Gray area mark the linear 
regime B < B* . 

the integration under these assumptions are depicted in 
Fig.[2j3, together with the subtracted diamagnetic contri- 
bution. 

Before discussing the results, let us first examine pos- 
sible sources of inaccuracy. The diamagnetic contribu- 
tion is small, and therefore can affect only the high- 
temperature susceptibility. The largest error may come 
from the extrapolation of dx/dn to zero density, as 
well from the uncertainty m 10 10 cm~ 2 in the density 
itself [ll|. We note that even if we make the extreme 
assumption that the magnetization is zero below til(T), 
thus significantly underestimating the susceptibility, we 
would get qualitatively similar results, as also shown in 
Fig.[2b. 

The most striking features of the susceptibility x are 
(i) its low-temperature value, at maxima exceeding the 
Pauli one by the factor of 40, and (ii) strong tempera- 
ture and relatively weak density dependence at high den- 
sities. The temperature dependence of the susceptibility 
per electron for several densities is shown in Fig.[3j It 
can be fit reasonably well with a power law x ^ 
with a w 2.4 ± 0.2. Not only does the susceptibility di- 
verge faster than the independent-spin Curie susceptibil- 
ity, n B /T, but its low-temperature low-density value ex- 
ceeds the Curie value, indicated by the solid line in Fig. [3] 
We emphasize that neither the susceptibility value, nor 
its temperature dependence change qualitatively when 
the system passes through the metal-insulator transition 
at n c ss 8.5- 10 10 cm -2 . Remarkably, the results for the R 
and H samples, shown in Figs. l2l3l bv the bold and hollow 
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FIG. 2: (a) Derivative of the low-field susceptibility with re- 
spect to the density dx/dn for different temperatures. Hollow 
symbols-R sample, bold symbols - H sample, (b) % obtained 
by integration; Dotted line (also in c) - non-renormalized 
Pauli susceptibility, dashed line-subtracted diamagnetic con- 
tribution, x - the same as □ integrated from the density 
til- (c) Susceptibility defined as rn/B at 2T. Solid line - the 
Shubnikov-de Haas susceptibility [SJ] . The temperatures are 
indicated by the symbols: □ - 1.7K; O -2.4K;< - 4K; O - 7K; 
* - f3K. Inset- maximum of the magnetization in the field 
< B < 2T at T = 1.7K, m and n are in units of 10 u cm~ a 

symbols respectively, are very similar. The largest con- 
tribution to the integral l|2|). which determines even at 
high densities comes from dx/dn at low densities. There- 
fore it is important to stress that dx/dn also diverges as 
1/T a , even deep in the metallic regime, as indicated in 
Fig.H 

The magnetic moment can be obtained similarly to 
the susceptibility, by integration of dm(n,T, B)/dn over 
n for given values of B and T. Since at low-temperature 
dm/dn saturates, and even drops for B > B*, as seen in 
Fig.DJt, m(B) behaves similarly; see the inset in Fig. 2] 
This saturation resolves the apparent contradiction be- 
tween the current thermodynamic and earlier 0,0] trans- 
port susceptibility measurements. Indeed, all the trans- 
port measurements to date were performed in the high- 
field domain, when the Zeeman energy is larger than 
the temperature, whereas we observe the divergent sus- 
ceptibility solely in the low field domain. As seen in 
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FIG. 3: Low-field susceptibility \ per electron vs. temper- 
ature for different densities, indicated in units of 10 n cm~ 2 . 
Hollow symbols-R sample, bold symbols - H sample. Also 
\dx/dn\ is shown for comparison with x and +. Dashed line 
line-x oc I/T 2 , solid line-Curie susceptibility. 

Fig. [4] B* depends linearly on temperature. Addition- 
ally, we note that B* appears to be density-independent 
for n iS 3 • 10 11 cm -2 . The low-field domain, at which 
the zero-field susceptibility can be determined, lies in 
the bottom-right corner. In contrast, both Shubnikov- 
de Haas, and magnetoresistance data belong to the top 
left corner. If one takes the magnetic moment at e.g. 2T, 
which is the typical total field in 0,0,0], ne would get the 
value of the susceptibility, defined as x(-B) = m(B)/B, 
similar to obtained in 0,0], as shown in Fig. [2b- 

In our previous work 0] we could only measure dm / dn 
down to « n c , and therefore we chose to integrate it 
from a high density at which we assumed m to be known 
and temperature independent. We used the susceptibility 
obtained from the Shubnikov-de Haas measurements in 
the high-field domain, (similar assumption was later used 
in [13[ ) as the initial value for the integration in both high 
and low field domains. As the present work shows the 
temperature dependence of the susceptibility in the low 
field domain cannot be ignored. 

A susceptibility diverging faster than 1/T indicates a 
ferromagnetic interaction between spins. If, however, the 
low temperature extrapolation of B* = ksT '/5.7 fig + 
0.03T in Fig. [4] is taken seriously, it means that the di- 
vergence should be cut off at about lOOmK. A divergent 
susceptibility contradicts numerical calculations [121 ] for 
a two-valley system, which predicts no tendency to fer- 
romagnetism. At high enough density one would expect 
the susceptibility to be determined by excitations in the 
vicinity of the Fermi level. Interaction corrections to sus- 
ceptibility in a clean system were calculated in 14 , [TH . 
Although the results differ in prefactors, both predict 
a correction to the susceptibility of the order of T/Ep- 
It is important to stress that not only docs x diverge 
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Magnetoresistance [4] 





Low-field domain: 
X~1/T a a=2.2-2.6 
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FIG. 4: Comparison with transport measurements. ■ marks 
almost density-independent B* . Arrows indicate the domains 
of the previous measurements of \- Inset: field dependence 
of magnetization for R sample at 1.7K. 

as 1/T a , but so does dx/dn, even deep in the metallic 
regime at densities n = 2 — 4 • 10 n cm~ 2 , in clear con- 



tradiction with [lj, [l_5j . We therefore think the divergent 
susceptibility is disorder-related and signals the presence 
of local moments in the 2D system up to high densi- 
ties; numerical calculations indeed predict susceptibility 
enhancement by disorder [l6[ ■ This suggestion does not 
contradict numerous existing transport data, since local- 
ized states are very rarely seen in transport, which probe 
the electrons on the time scale of ps. In contrast, ther- 
modynamic measurements probe all the states which can 
be recharged on the time scale set by the magnetic field 
modulation. On the other hand, it is important to em- 
phasize that the maximal magnetization of the system, 
shown in the inset in Fig.[2p,, is very close to n/xs at low 
densities. This rules out the contribution of distant mag- 
netic centers, existing in addition to the 2DEG (such as, 
e.g., localized states in the bulk Si, or deep levels at the 
interface). 

A clue about the origin of the anomalous susceptibility 
may come from the magnetization drop at low temper- 
atures in the intermediate magnetic field, above B* and 
below some 2.5T, visible in the inset in Fig.[4j the magne- 
tization grows again at even higher fields. For a magnetic 
field coupled only to the spins, the magnetization should 
grow with the field. Therefore this drop must be related 
either to the orbital effects, or to the field induced varia- 
tion in the number of electrons contributing to the signal. 

In conclusion, we observed a strongly enhanced diver- 
gent susceptibility at low temperatures in high-mobility 
MOSFET's. We find no qualitative change of the sus- 



ceptibility behavior in the vicinity of the metal-insulator 
transition density n c , and the divergence persists deep 
into the high density metallic phase, similar to the ob- 
servations on phosphorus doped Si. 
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Programs of RAS, Russian Ministry for Education and 
Science, and the Program "Leading scientific schools". 
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